Origin of Hysteresis in a Proximity Josephson Junction 
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We investigate hysteresis in the transport properties of Superconductor - Normal metal - Super- 
conductor (S-N-S) junctions at low temperatures by measuring directly the electron temperature 
in the normal metal. Our results demonstrate unambiguously that the hysteresis results from an 
increase of the normal metal electron temperature once the junction switches to the resistive state. 
In our geometry, the electron temperature increase is governed by the thermal resistance of the 
superconducting electrodes of the junction. 

PACS numbers: 74.50.+r, 74.45,+c 



Dissipationless supercurrent can flow between two su- 
perconductors, up to a critical current I c in a Joseph- 
son junction. The junction dynamics can be described 
by the Resistively and Capacitively Shunted Junction 
(RCSJ) model [1]. The junction capacitance is then re- 
sponsible for hysteresis in the current- voltage character- 
istic. In lateral junctions, the distance between the two 
superconducting electrodes induces an extremely small 
capacitance, much lower than those in a typical tunnel 
junction. A non-hysteretic (overdamped) current- volt age 
characteristic is then expected. Nevertheless, a signifi- 
cant hysteresis is routinely observed in lateral junctions 
as soon as their critical current is large: once the junc- 
tion has switched to the resistive branch, it does not re- 
cover the superconducting state until the bias current 
is decreased to a significantly smaller retrapping current 
I r . This observation does not depend on the nature of 
the weak link, as it was early observed in superconduct- 
ing constrictions and microbridges 0, EjjL 0], an d more 
recently in superconducting nanowires [5], normal met- 
als 0,13? two-dimensional electron gases (2-DEG) Q, 
semiconductor nanowires [9], carbon nanotubes [Io[ and 
graphene [11]. Two main explanations have been pro- 
posed. Firstly, the Joule power deposited in the weak 
link can induce a self-heating process so that the local 
temperature in the normal part increases [2|, y, ll2[ . Sec- 
ondly, it has been proposed that the response time of the 
junction R n C should be replaced by a time h/A related 
to the superconducting gap A [4] or by the electron diffu- 
sion time L 2 /D through the junction @. Here R n is the 
normal state resistance, L is the junction length and D is 
the electron diffusion constant. These latter approaches 
are equivalent to assuming an effective capacitance larger 
than the geometric capacitance. Despite relatively good 
agreement with the experiments, no strong justification 
can be brought to support these hypotheses. 

In this Letter, we report an experimental study of 
the hysteretic regime of proximity Superconductor - Nor- 
mal metal - Superconductor (S-N-S) Josephson junctions. 
Our results demonstrate unambiguously that the hystere- 
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FIG. 1: (Color online) Top: SEM image of sample 2 contain- 
ing a S-N-S junction of 1.5 /mi length with a sketch of the 
measurement circuit. Two tunnel probes (top of the image) 
are connected to the normal metal embedded between two 
superconducting banks (on left and right sides of the image). 
The overlap of the superconducting banks (dark grey) with 
the normal metal layer (light grey) is visible. During the mea- 
surement, the S-I-N-I-S junction is biased at a fixed current 
Isinis and the voltage drop Vsinis is monitored. Bottom: 
Temperature dependence of the voltage Vsinis of the same 
sample with a current bias Isinis = 6 pA (no current is 
flowing through the S-N-S junction)and zoom of the sample 
image. 



sis results from the increase of the normal metal electron 
temperature once the junction switches to the resistive 
state. An electron temperature of up to 0.6 K is mea- 
sured while the thermal bath remains at 50 mK. We show 
that, in our geometry, the electron temperature increase 
is governed by the thermal resistance of the supercon- 
ducting electrodes. 

Figure [1] shows a SEM image of one of the three sam- 
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# L (/im) i2n (0) 


w (nm) 


E T h (/ieV) Ethju (^eV) 




0.0 


1 1.0 6.32 


100 


4.8 3.9 


0.65 




2 1.5 9.96 


100 


2.4 2.1 


0.48 


0.1 


3 2.0 8.46 


150 


1.4 2.6 


0.040 





TABLE I: Sample parameters. The normal state resistance 
R n , the distance between the two superconducting electrodes 
of the S-N-S junction L, the N wire width w and the esti- 
mated Thouless energy E T h taking into account the junction 
length plus the overlaps of 0.1 /im with the superconduct- 
ing electrodes are listed. Fit parameters of Fig. [3] are the 
Thouless energy Ethju and the reduction parameter a. The 
superconducting gap A is about 200 fieV for each sample. 



pies we have measured. It consists of a S-N-S junction 
with two additional superconducting tunnel probes on 
the normal metal part of the junction, thus forming a S- 
I-N-I-S junction (I stands for Insulator) sharing its nor- 
mal metal part with the S-N-S junction. The sample was 
made by three-angle deposition through a suspended re- 
sist mask. Every structure is then tripled, so that the 
Josephson junction leads are actually made of the over- 
lap of two superconducting (Al) layers and one normal 
metal (Cu) layer. The two tunnel probes were realized 
by depositing Al as the first layer and oxidizing it in 
a pressure of 5 mbar for 5 minutes prior to Cu depo- 
sition. The tunnel junction resistances are then in the 
range of 30-100 kQ each. A thick (70 nm) Al deposition 
was finally performed in order to obtain the supercon- 
ducting electrodes of the S-N-S junction. Based on our 
measurements, the residual resistance of the N-S contacts 
is estimated to be below 0.5 fi. The parameters of the 
samples are listed in Table 1 . From the measured normal- 
state resistances, the diffusion constant of Cu is inferred 
to be about 110 cm 2 /s for every sample. The normal 
metal (Cu) strip width and thickness are about 100-150 
nm and 27 nm respectively. The junction length is 1, 
1.5 and 2 /im for sample 1, 2 and 3 respectively. Here, 
the junction length L is much larger than the supercon- 
ducting coherence length £ s , so that the Thouless energy 
Eth — TiD / L 2 [14[ is much smaller than the supercon- 
ducting gap A. In this so-called long junction limit, the 
Thouless energy defines the magnitude of the proxim- 
ity effects 14 , including the proximity-induced mini-gap 
width 15|, LL6[ and the critical current [17[ . 

Tunneling through a S-I-N-I-S junction is sensitive to 
the electron temperature in the normal metal [l3j |. Here, 
we bias the S-I-N-I-S junction with a battery-powered 
current source in the pA range and measure the volt- 
age, see Fig. [T] top. Figure [T] bottom left shows the 
temperature dependence of the sample 1 S-I-N-I-S volt- 
age at a fixed bias current of 6 pA, with zero current 
applied through the S-N-S junction. We obtain the ex- 
pected almost linear behavior in a wide temperature 
range, without saturation down to below 50 mK. With 
this calibration, the S-I-N-I-S junction serves as an elec- 
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FIG. 2: (Color online) Current- volt age characteristic of sam- 
ple 1 S-N-S junction (bottom) shown on the same current 
scale with the S-I-N-I-S thermometer voltage response (top) 
measured simultaneously at a 50 mK cryostat temperature. 
In the top part, the right vertical axis gives the corresponding 
electron temperature. 



tron thermometer. In the following, we will assume an 
electron population in the normal metal close to quasi- 
equilibrium, so that we can define an effective electron 
temperature T e , possibly different from the cryostat tem- 
perature T ba th- 

Figure [2j bottom part, displays the S-N-S junction 
current- voltage characteristic of sample 1 at the cryostat 
base temperature (Tb at h = 50 mK). It features a clear 
superconducting branch at zero voltage. When the cur- 
rent is increased, a sudden switch to a resistive branch 
with a constant resistance can be seen. As previously dis- 
cussed, the characteristic is hysteretic. When the current 
is decreased, the voltage jumps back to zero only at a re- 
trapping current significantly smaller than the switching 
current. Like the switching, the retrapping appears as 
a discontinuity of the characteristic. Figure [2] top part 
displays the simultaneously measured voltage response 
of the electron thermometer at a 20 pA current bias. 
Here we have subtracted the voltage drop in the nor- 
mal metal, between the two tunnel junctions, due to the 
S-N-S bias current. The corresponding electron temper- 
ature T e scale is given. A striking behavior is observed. 
In the supercurrent branch, the electron temperature is 
almost constant as expected. Yet it changes slightly be- 
cause of heating by the current through resistive filters on 
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FIG. 3: (Color online) Temperature dependence of the R n Ic 
product, where I c is the measured critical current, for samples 
1-3 (open symbols). The fits are displayed as full lines. For 
sample 3, the theoretical prediction using a = 0.5 and Eth — 
1.12 /xeV (close to the expected value) is shown as a dotted 
line. The R n I r product including the retrapping current I r 
at a bath temperature of 50 mK is plotted for each sample 
versus the electron temperature just before retrapping (full 
symbols, indicated by arrows). 



the sample stage. At the switch to the resistive branch, 
the electron temperature jumps to a much higher value. 
After the jump, the electron temperature still increases 
because of the increased Joule power. When the current 
is decreased, the electron thermometer signal first follows 
the same curve. Temperature stays elevated until it drops 
to the bath temperature value, precisely at the retrapping 
current. This demonstrates clearly that the hysteresis in 
our S-N-S junctions is governed by over-heating of the 
normal metal. The same measurement was performed on 
the same sample 1 at a different current bias of 12 pA, 
giving the same electronic temperature evolution, and on 
samples 2 and 3, displaying similar behavior. 

Before turning to a quantitative analysis of our experi- 
mental data, let us discuss the possible out-of-equilibrium 
effects. This will be done by comparing the retrapping 
and the switching current data. The measured temper- 
ature dependence of the R n I c product for each sample 
is plotted in Fig. [3l We fitted the data to the theo- 
retical prediction [17[ using as fit parameters the Thou- 
less energy and a scaling parameter a accounting for the 
non-ideality of the N-S interfaces. For the two shorter 
samples (1, 2), we obtain a good fit with values of the 
Thouless energy close to the estimates and values for a of 
about 0.5 [18]. Although fabricated on the same chip as 
sample 2, sample 3 shows a different behavior with a re- 
duced critical current and an increased effective Thouless 
energy, which we interpret as due to the finite phase co- 
herence length compared to the junction length. In Fig. 
[3j we also plot the R n I r product versus the electronic 
temperature before retrapping (square symbols). Here 
I r is the retrapping current at a bath temperature of 50 



mK. In a quasi-equilibrium hypothesis, retrapping occurs 
when the bias current is equal to the critical current at 
the electron temperature so that the latter points should 
coincide with the equilibrium critical current data. The 
retrapping data are actually close to the equilibrium data 
with a shift of about 50 mK or less than 50 % in critical 
current amplitude. This limited discrepancy shows that 
a discussion in terms of effective electron temperature is 
reasonably justified. 

An elementary idea to analyze the electron heating is 
to consider the hot electrons to be confined in the normal 
metal. In this case, the electron-phonon interaction en- 
sures the coupling to the thermal bath pjl • In the limit of 
a low temperature for the phonons, the power flow writes 
P = £/7T e 5 with £ = 2 nW./im- 3 .K- 5 in Cu and U is 
the metal volume. Here, this power flow is equal to the 
Joule power IsnsVsns dissipated in the S-N-S junction. 
In the case of sample 2 at 1 fiA injected current, the pre- 
dicted quasi-equilibrium temperature T e = (P/EZ7) 1 / 5 is 
about 1 K. Figure [4] displays, as dotted lines, the power 
calculated within this hypothesis for the three different 
samples, in parallel with the experimental data, as a func- 
tion of the inverse of the measured electronic tempera- 
ture T e . At a given temperature, the experimental data 
is well above the prediction, which means that electrons 
thermalize via another process. 

We consider now the thermal link through the super- 
conducting electrodes attached to the normal metal. The 
two tunnel N-I-S junctions are expected to be good ther- 
mal insulators compared to the transparent N-S inter- 
faces. In the superconducting state, quasi-particles with 
an energy above the superconducting gap contribute to 
the heat transport through the interface and the super- 
conductor. This was recently discussed in the framework 
of noise measurements in a S-N-S junction in the hot elec- 
tron regime [20[. In the following, we neglect the thermal 
resistance of the S-N interface compared to that of the 
superconducting line. The ratio r between the thermal 
conductivities in the superconducting state and in the 
normal state writes [21[ 

The normal-state thermal conductivity is = Log^T, 
with gN being the normal state electrical conductivity. 
Here, we assume the Wiedemann- Franz law with the 
Lorentz number L =2.45.10 -8 W.^.K" 2 . In the low 
temperature limit, the quantity r(T) decays exponen- 
tially with the temperature as exp(— A/kBT e ). 

The power flow as a function of the electronic tem- 
perature can then be calculated as an integral of the 
temperature-dependent thermal conductivity: 

P(T e ) = L G% / r(T)TdT (2) 

JT hath 
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the dissipated power density at the switching is larger. In 
the case of Josephson junctions based on nanowiresJHI, 0] , 
2-DEGs carbon nanotubes or graphene [ll|7the 
power density at the switching is estimated to be about 
or above 1 nW//im 3 , while the samples investigated here 
feature generally a smaller density in the range 2.10- 3 to 
1 nW//im 3 . This suggests that the observed hysteresis 
in these other kinds of lateral Josephson junctions is also 
due to electron heating. 
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FIG. 4: (Color online) Injected Joule power IsnsVsns in the 
S-N-S junction as a function of the inverse of the measured 
electronic temperature T e . The cryostat temperature is 50 
mK. The expected behavior due to the electron-phonon cou- 
pling is plotted as dotted lines. The dash-dotted line is the 
power calculated using Eq. [2] and A = 200 fieV. The numbers 
refer to the samples. 



where is the normal state electrical conductance of 
the superconductor from the hot region to the thermal 
bath. Here, we expect the Cu film to be well thermalized 
and act as a thermal bath. The relevant conductance is 
thus determined by the 70 nm-thick Al electrode between 
the N-S interface and the overlap with the Cu film. This 
corresponds to a resistance of about 9 square resistances 
on both electrodes (see Fig. [1]). Assuming a typical re- 
sistivity of 2 /ifi.cm for Al, this gives G s ^ ~ 1.3 Vt. The 
related normal-state thermal conductance LqG s ^ is then 
about 19 nW/K at 1 K. 

The power calculated with this parameter is displayed 
as a function of the inverse of the electronic tempera- 
ture T e in Fig. 2] (dash-dotted line). We took for the 
S-N-S junction a superconducting gap value A = 200 
fieV equal to the one measured in the S-I-N-I-S junction. 
The calculation compares favorably with the experimen- 
tal data. On a semi- logarithmic plot, both show a nearly 
linear decay of the power with the temperature inverse. 
In most of the investigated range, the measured power is 
lower than the prediction, which means that the thermal 
conductance of the superconductor was over-estimated. 
Nevertheless, the fair agreement between the data and 
our simple model shows that the electron temperature 
increase is actually limited by the thermal conductance 
of the S-N-S junction superconducting electrodes. 

In conclusion, our study solves a long-standing issue in 
the general field of Josephson junctions by showing that 
the hysteresis routinely observed in long S-N-S junctions 
is of thermal origin. This interpretation should definitely 
hold also in shorter metal-based S-N-S junctions where 
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